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^ We have studied the infrared properties of the brightest cluster galaxies 

(BCGs) located in the cores of X-ray-luminous clusters at 0.15 < z < 0.35. 
^ ! The majority of the BCGs are not particularly infrared-luminous compared with 

other massive early-type galaxies, suggesting that the cluster environment has lit- 
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tie influence on the infrared luminosities of the BCGs. The exceptions, however. 



', are the BCGs in the three X-ray-brightest clusters in the sample, A1835, Z3146, 

5^ I and A2390. These BCGs have a prominent far-infrared peak in their spectral 

^ I energy distributions (SEDs), and two of them (those in A1835 and Z3146) can 

\^ I be classified as luminous infrared galaxies (LIRGs: Ljji > 10^^ L©). Although 



radio AGNs are found to be prevalent among the BCGs, the infrared luminosi- 
ties of these three BCGs, judged from the infrared SED signatures, are likely to 
be powered by star formation. Considering the overall trend that clusters with 
^ \ shorter radiative gas cooling times harbor more infrared-luminous BCGs, the en- 

^ I hanced star formation may be caused by the cooling cluster gas accreting onto 

the BCGs. 
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Subject headings: galaxies: clusters: general — cooling flows — galaxies: cD — 
galaxies: active — infrared: galaxies 

1. Introduction 

Clusters of galaxies are embedded in a halo of diffuse hot gas {Tgas ~ 10^ — 10^ K), which 
emits strongly in the X-ray through thermal bremsstrahlung radiation. The radiative cooling 
time of this hot intracluster plasma is usually comparable to the Hubble time, suggesting that 

the gas does not cool substantially over the age of the clusters. There are, however, clusters 
that show a highly peaked X-ray surface brightness profile in the core, which indicates that 
the hot gas in these cores is dense and has a substantially shorter cooling time. In extreme 
cases, the cooling time is estimated to be less than a few Gyrs, which prompted the idea that 
these dense cores are sustaining an inflow of cooling gas toward the cluster center (Cowie & 
Binney 1977; Fabian & Nulsen 1977). 

The reality of such a cluster "cooling flow" has been disputed over the years on the 
basis that the inferred gas cooling rate (which we assume to lead to a mass deposition rate) 
is simply too large to be consistent with the modest level of activity (e.g., Ha-derived star 
formation rates) seen in these cluster cores (see Fabian (1994) for a review of cooling flows). 
In the most extreme cases, the gas coohng rate was estimated to be >1000 M0 yr~^ (e.g., 
Allen et al. 1996), and yet there is no sign of activity involving such a magnitude of mass 
inflow. Many of the brightest cluster galaxies (BCGs) in strongly cooling cores do show some 
signs of activity, such as visual emission lines and UV continuum (e.g., Allen 1995) , but 
if they are due to star formation, the inferred star formation rate (e.g., from Ha) is only a 
few hundred Mq yr~^ at most (e.g., Crawford et al. 1999). Considering that the BCGs in 
strongly cooling cores also show evidence for containing a signiflcant amount of dust (Edge 
et al. 1999) and molecular gas (Edge 2001; Salome & Combes 2003), it is clear that these 
BCGs are unusually active; what is puzzling, however, is the large difference between the 
gas cooling rate we infer and that we can account for. 

A self-consistent picture is flnally emerging from observations with Chandra and XMM. 
XMM especially has demonstrated convincingly for the first time that the cluster gas is 
coohng but not as much as originally thought. More specifically, XMM has detected strong 
emission from a cool plasma at temperatures half to one-third of the ambient value, but 
there is also a severe deficit of emission from plasma below these temperatures (Peterson 
et al. 2001; Tamura et al. 2001; Kaastra et al. 2001; Peterson et al. 2003). In other words, 
there seems to be a floor of gas temperature below which cooling is no longer effective. The 
result is an order of magnitude reduction in the X-ray derived mass deposition rates, making 
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them consistent with the modest amount of activity observed at other wavelengths, such as 
the star formation rates derived from the extinction-corrected Ha emission (Crawford et al. 
1999). 

The current view is that some kind of gas heating (e.g., a radio jet from an AGN) slows 
the coohng process (see Begelman (2004) for a recent review). Indeed, such heating may 
play an important role in the process of galaxy formation in general because it would provide 

a way to suppress the formation of bright (> L*) massive galaxies, which are overproduced 
in model calculations (Croton ct al. (2006) and references therein). In fact, the BCGs in the 
cooling-fiow clusters may prove to be an important laboratory where we can study the effects 
of cooling gas directly accreting onto a seed mass concentration, a condition that might be 
found in the first generation of forming galaxies at high redshift. 

To investigate this cooling flow problem further, we have examined the Spitzer/YRAC 
(3.6, 4.5, 5.8, and 8.0 /xm) and Spitzer/MlV^ (24, 70, and 160 /xm) images of the BCGs 
located in the cores of X-ray-luminous clusters. At the redshifts of targeted clusters, IRAC 
detects stellar continuum hght while MIPS detects thermal dust emission. Observations of 
emission from heated dust can detect deeply embedded forms of activity that are hidden 
from other search methods such as visual emission lines or UV continuum. For example, a 
galaxy with a star formation rate of ~ 1000 yr~^ would be ultraluminous in the infrared 
(LiR > 10^^ L0), so if any of the BCGs in the so-called cooling flow clusters are undergoing 
such vigorous star formation as a result of a massive gas inflow, they would stand out in the 
mid-/far-infrared images. The overall shape of spectral energy distribution (SED) 
will also tell us whether the infrared luminosity source is star formation or an AGN. 

Throughout the paper, we adopt the cosmological parameters of VLm — 0.3, Oa — 0.7, 
and Hq = 70 km s~^ Mpc~^ except when noted otherwise. 

2. The Sample 

The eleven clusters for this study are listed in Table 1 and were part of a Spitzer GTO 
program described in Egami et al. (2006). The overall cluster selection criteria for this 
program are the following: (1) X-ray luminous {Lx ^ 5 x 10'^^ erg/s based on the published 
soft X-ray luminosities listed in Table 1), (2) moderate redshift {z ~ 0.15 — 0.5), and (3) 
low IR background {Nh < 3.5 x 10^^ cm~^). We did not include clusters whose redshifts are 
based on a small number of measurements and those with severely limited ancillary data. 
Since the primary goal of this program is to study gravitationally lensed background galaxies, 
the most X-ray-luminous clusters were chosen under the assumption that they are the most 
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massive clusters and therefore the most effective lenses. This means that these clusters were 
chosen without any regard to the properties of the BCGs. However, since, to a first order, 
the X-ray luminosity is proportional to the gas cooling rate, this cluster sample also provides 
a good data set to study the properties of BCGs in strongly cooling cluster cores. 

The identification of the BCG in a cluster is usually straightforward. Many cluster 
cores harbor a single giant elliptical galaxy that dominates the brightness at the visual/near- 
infrared wavelengths and is located at/ncar the peak of the cluster X-ray emission. There 
are, however, cases in which a BCG cannot be identified uniquely: some clusters have more 
than one giant elliptical with similar brightnesses in the chister core, or have a disturbed 
X-ray morphology without a well-defined peak. In these cases, it was not clear how we could 
relate the cluster gas properties (e.g., cooling rates) with those of the individual galaxies. 
For this study, we therefore chose only those clusters that show an X-ray morphology with a 
well-defined single peak and have a single 24 //m-detected BCG near the X-ray peak based 
on the archived ROSAT HRI images and on published Chandra images (e.g., Smith et al. 
2005). We also note that based on the C/ian(/ra images, our sample contains both "relaxed" 
and "unrelaxed" clusters, and that the spatial coincidence of the X-ray peak and BCG is 
not precise in many of the "unrelaxed" clusters (e.g.. Smith et al. 2005). 

Table 1 lists the position of the BCG in each cluster as well as a few critical properties 
of these clusters such as X-ray luminosities and gas cooling times. Based on the classification 
by Allen & Fabian (1998), 7 out of the 11 clusters are classified as cooling flow clusters (CF) 
while four are classified as non-cooling fiow clusters (NCF). 

3. The Data 

3.1. Spitzer/IRAC Data 

The 3.6, 4.5, 5.8, and 8.0 /xm images were obtained using the Infrared Array Camera 
(IRAC; Fazio et al. 2004) on board the Spitzer Space Telescope (Werner et al. 2004). Each 
of the four IRAC channels uses a separate detector array. The 3.6 iim channel (Ac = 3.56//m; 
AA = 0.75/xm) and 4.5 iim channel (Ac = 4.52//m; AA = 1.01//m) use 256 x 256 InSb arrays 
while the 5.8 /im channel (Ac = 5.73//m; AA = 1.42//m) and 8.0 /im channel (Ac = 7.91/im; 
AA = 2.93//m) use 256 x 256 Si:As arrays. At any given time, two channels are used 
simultaneously (3.6 and 5.8 fira; 4.5 and 8.0 /im). Each channel has a similar pixel size 
(~1'.'2 pixel"^) and field of view (5f2x5f2). Integration times were 2400 s per band except 
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for A1835, whose integration times were 2400 s at 3.6 & 5.8 /xm and 3600 s at 4.5 & 8.0 //m^. 
The images were taken with the small-step cycling dither pattern. The basic calibrated data 
(BCD) images were combined using a custom IDL mosaicking routine with the final pixel 
scale of 0'.'6 pixel"^, half of the instrument intrinsic pixel size. 

The IRAC photometry used a circular beam with a diameter of 12'.'2 with a sky back- 
ground annulus of 6'.'1-12'.'2 in radius. We selected this rather large beam (6x the PSF 
FWHM) because at these wavelengths, some BCGs show a significant spatial extent. Al- 
though in some cases there is contamination from nearby sources falling into the beam, the 
effect is small because of the dominance of the BCGs. The point-source aperture correc- 
tions were applied, which were 1.061, 1.064, 1.067, and 1.089 at 3.6, 4.5, 5.8, and 8.0 /xm, 
respectively, based on the IRAC Data Handbook. 

The resultant 1 a sensitivities in the mosaicked images are ~ 5/iJy in all four bands. 
The shorter-wavelength images do not necessarily go deeper when such a large photometry 
beam is used with IRAC images of this depth due to the confusion noise. A conservative 
estimate for the absolute calibration uncertainty is 10%. 

3.2. Spitzer/MlVS Data 

The 24, 70, and 160 ^m. images were obtained using the Multi-band Imaging Photometer 
for Spitzer (MIPS; Rieke et al. 2004). The MIPS 24/im channel (Ac = 23.7/xm; AA = 4.7/xm) 
uses a 128x128 BIB Si: As array with a pixel scale of 2'/55 pixel"^. The 70 jim channel 
(Ac = 71.4//m; AA = 19.0//m) uses a 16x32 Ge:Ga array with a pixel scale of 9'.'98 pixel"^. 
The 160 jjim channel (Ac = 155.9/im; AA = 35//m) uses a 2x20 stressed Ge:Ga array with a 
pixel scale of 16"xl8". 

All the clusters in the sample (Table 1) were imaged at 24 //m using the large-source 
photometry mode. The total integration times were 3600 s pixel"^ at the position of the 
BCGs. Since we did not have enough observing time to observe all the BCGs at 70 and 160 
//m, at these wavelenghts we only observed the three BCGs brightest at 24 /xm, for which 
we can expect a good chance of detection. For the two clusters with the brightest BCGs 
at 24 yum (A1835 and Z3146), we conducted 70 and 160 ^va. imaging using the large-source 
photometry mode. For the cluster with the third brightest BCG at 24 yum (A2390), we 
performed only 70 yum imaging due to a high infrared sky background toward this direction. 



iThe IRAC data of A1835 were taken as part of the GTO program PID:64 (PI: Giovanni Fazio). All the 
other data were taken as part of the GTO program PID:83 (PI: George Rieke). 
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The integration times were 370 s pixel ^ at 70 /xm, and 100 s pixel ^ at 160 /xm at the 
position of the BCGs. 

The data were reduced and combined with the Data Analysis Tool (DAT) developed 
by the MIPS instrument team (Gordon et al. 2005). A few additional processing steps were 
applied to the 24 //m data as well. First, special fiat fields were constructed to remove 
artifacts in the flats produced by material deposited on the scan mirror. Secondly, the 
24 ^m array is subject to latent images after exposure to a bright source. The latent images 
manifest themselves as a region of reduced sensitivity centered on the array position of the 
bright source in subsequent exposures. These "dark latents" decay with time and disappear 
upon the next anneal cycle. Several of our observations were contaminated by dark latents 
caused by bright 24 /xm sources in the target fields. A source-free image of the latents was 
produced by median combining the entire data set with appropriate bright source masks. 
The dark latents were removed from each individual frame by dividing the frames by the 
normalized dark latent image prior to mosaicking. A detailed discussion of these instrumental 
effects will be presented in a forthcoming paper (G. H. Rieke et al. 2006, in preparation). A 
third processing step beyond the normal DAT processing was necessary owing to the design 
of the observations. Each cluster was observed using two Astronomical Observing Requests 
(AORs), and the two AORs were often separated by a significant amount of time. As a 
results, changing zodiacal backgrounds resulted in strong gradients in the resulting mosaics. 
To mitigate the effect of any DC offset between AORs, a mean row was subtracted from each 
individual data frame, again with an appropriate bright source mask, prior to mosaicking. 
The final mosaicked images were produced with pixel scales of l'/25, 4'/92, and 8" pixel"^ at 
24, 70, and 160 /im, respectively, which are half of the instrument intrinsic pixel sizes. 

The MIPS photometry beam diameters and corresponding inner/outer sky annulus radii 
are as follows: 26" and 20"-32" at 24 //m, 35" and 39"-65" at 70 //m, and 50" and 75"-125" 
at 160 /im. The point-source aperture corrections for these measurements are 1.167, 1.308, 
and 1.442, respectively, based on the MIPS instrument web page. For the clusters A2218 and 
AC114, we used a smaller photometry beam (13'.'2 diameter, aperture correction = 1.648) 
at 24 yum to prevent contamination from nearby sources. The 70 and 160 fim aperture 
corrections depend to some extent on the intrinsic color of the observed source, so we used 
the value for a red PSF (T = 15 K) at 70 /xm and that for a blue PSF (T = 1000 K) 
at 160 iim, assuming that the infrared SEDs of the galaxies observed here are peaking 
somewhere between these two wavelengths. However, the effect of this color-dependent 
aperture correction is on the order of a few to several percent at most, well within the 
uncertainties of our photometric measurements. 

The resultant 1 a sensitivities in the mosaicked images are 0.01, 2.6, and 15 mJy at 24, 
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70, and 160 //m, respectively. The absolute calibration uncertainties are 10, 20, and 20% in 
these three bands. 



3.3. Chandra X-ray Data 

For the three X-ray brightest clusters in the sample (A1835, A2390, and Z3146), we have 
extracted all available data from the Chandra archive, and have reprocessed them to provide 
a uniform and consistent set of X-ray measurements. The initial Chandra results for A1835 
have been published previously in Schmidt, Allen, & Fabian (2001). For A2390, an analysis 
of the combined 19 ksec observation from Cycle 1 has been published in Allen, Ettori, & 
Fabian (2001); however, a subsequent deeper exposure was taken in Cycle 4, and we also 
included this data set. Finally, although observed during Cycle 1, the Chandra results for 
Z3146 remain unpublished. When all the data were combined, the total integration times for 
A1835, A2390, and Z3146 (after standard filtering) were 30 ksec (with two data sets), 104 
ksec (with three data sets) , and 46 ksec (with one data set) , respectively. The data processing 
and the method for calculating the gas coohng rates are described in Appendix A. Detailed 
analyses of the A2390 and ZW3146 Chandra data will be published elsewhere (M. W. Wise 
2006, in preparation). 

3.4. VLA 1.4 GHz Data 

We have compiled the 1.4 GHz flux densities of the BCGs from a number of sources. 
Where duplications exist, we have adopted the following priorities: 

1. Deep high- resolution observations targeting clusters (Morrison et al. 2003; Rizza et al. 
2003) 

2. FIRST survey (2003 Apr 11 version) (Becker, White, & Helfand 1995; White et al. 
1997) 

3. NVSS survey (Condon et al. 1998) 

The data published by Morrison et al. (2003) and Rizza et al. (2003) have the highest 
spatial resolution (I'.'S, A configuration) and sensitivity (cr ~ 25 — 50/xJy). The FIRST and 
NVSS surveys have spatial resolutions of 5" (B configuration) and 45" (D configuration), 
and sensitivities of u ~ 140 and 450 //Jy, respectively. The fiux densities from the NVSS 
should be used with caution because we see cases in which a single NVSS source at a BCG 
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position breaks up into multiple radio components (e.g., Z3146) or is clearly displaced from 
the BCG (e.g., A2219) when observed with higher spatial resolution. 

All the VLA flux densities are hsted in Table 2. When there is no radio detection (i.e., 
A773, A2218, and A2219), we list a 3 (t upper limit. 

4. Results 

Figure 1 and 2 show the IRAC and MIPS images of the core regions of the eleven clusters. 
The measured Spitzer&ux. densities are given in Table 2. Figure 1 shows the <Spzfeer images 
of the three BCGs that are the brightest at 24 //m. Figure 2 shows the rest of the sample. 
Note that these three BCGs are also the brightest in terms of the intrinsic mid-infrared 
luminosity since all the clusters studied here are at similar redshifts (2; ~ 0.2 — 0.3). The 
SEDs produced from the IRAC/MIPS photometry are shown in Figure 3^. 

The strong infrared sources in these three BCGs are compact, the 24 //m sources es- 
sentially being unresolved, which sets an upper limit on the source diameter (D) of D < 20 
kpc for the size of the infrared emitting region. The BCGs are resolved at 8 //m with an 
estimated intrinsic size of I? ~ 8 kpc, but the stellar light contribution (as opposed to ther- 
mal dust emission) may be significant at this wavelength. Therefore, we treat this size as an 
upper limit. This indicates that the infrared luminosity originates in the BCGs. 

The rest of the BCGs in the sample show SEDs monotonically decreasing in the mid- 
infrared with no clear indication of having a significant amount of far-infrared luminosity. 
These SEDs could still peak in the far-infrared, but the corresponding far-infrared lumi- 
nosities would be much smaller than those of the three infrared-bright BCGs. At the same 
time, all the BCGs are detected above 5 ex at 24 /im, indicating that these BCGs are not 
completely infrared-quiet. 

Figure 3 also shows that none of the BCGs has a power-law SED in the infrared typical 
of an AGN-dominated galaxy. Instead, all of them show a gradual drop at Ajcst 3 — 5fim, 
a characteristic of star-dominated SED (here, we see only the longer wavelength drop from 
the SED peak near 1.6 /im in the restframe). In fact, despite the large variation in the 
mid-/far- infrared range, the SEDs of the BCGs are quite similar at Arest ~ 3 — 5//m not only 
in terms of shape but also of luminosity. This indicates that the near-infrared light of these 
BCGs is dominated by equally massive stellar populations, which are not directly related to 



^In this and subsequent figures in the paper, error bars are omitted when their sizes are comparable to 
or smaller than those of the plot symbols. 
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the generation of the infrared luminosity. 

5. Discussion 

5.1. Infrcired Properties of the BCGs 

Figure 4 plots the restframe 15 /im monochromatic luminosities of the BCGs against 
the soft X-ray (0.1-2.4 keV) luminosities of their parent clusters. The soft X-ray luminosity 
is a good measure of the intracluster gas cooling rate, and unlike the hard X-ray, it will not 
be affected severely by AGN emission. Tabic 1 lists the i? 05/1 T soft X-ray (0.1-2.4 keV) 
luminosities from two i? T cluster catalogs, the i?(95'^T Brightest Cluster Sample (BCS) 
(Ebeling et al. 1998) and the Northern ROSAT All-Sky (NORAS) Galaxy Cluster Survey 
(Bohringer et al. 2000). As seen in Table 1, the luminosities listed in these two catalogs are 
consistent with each other with a ~10% uncertainty; we used the latter catalog by default. 
The two southern clusters, MS2137.3-2353 and AC114, are not listed in the ROSAT-ESO 
Flux Limited X-ray (REFLEX) Galaxy Cluster Survey (Bohringer et al. 2004), the southern- 
sky version of NORAS, so instead we used the Einstein 0.3-3.5 keV (Gioia & Luppino 1994) 
and ASCA 2-10 keV (Allen 2000) X-ray luminosities, respectively, as surrogates for the 
ROSAT soft X-ray luminosities. The latter needs to be divided by a factor of 1.25 to 
account for the systematic offset between the ASCA and ROSAT luminosities^. 

We have estimated the restframe 15 fim monochromatic luminosities based on the SEDs 
shown in Figure 3. We chose to calculate the monochromatic luminosity at 15 fim because 
the conversion from the 15 fim luminosity to the total infrared luminosity is well established 
based on the ISO observations (Elbaz et al. 2002). For comparison, the 1 a range of the 
15 fjim luminosities measured for local early-type galaxies by ISO (Ferrari et al. 2002) is 
shown (the two horizontal dotted hues). The figure shows that the majority of the BCGs 
studied here are not particularly luminous in the mid-infrared compared with the sample of 
Ferrari et al. (2002). Ferrari et al.'s sample may be biased toward infrared- active early-type 
galaxies, but this comparison at least shows that there is nothing special about the infrared 
luminosities of most BCGs, suggesting that the cluster environment has little influence on 
the BCG infrared properties in general. 



^Although the ASCA passband covers the hard X-ray and not soft X-ray, Table 1 shows that emperically 
there is a good scaling relation between the ASCA and i?05'j4T luminosities in the sense that the former is 
1.25 ± 0.05 times larger than the latter when A2219 and A2390 are exluded. The hard X-ray luminosities of 

these two clusters arc likely boosted by AGN emission. As we will show later, the BCG in A2390 contains 
a strong radio AGN while A2219 is known to have a strong radio source in the cluster core near the BCG. 
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The only exceptions are the BCGs in A1835, Z3146, and possibly A2390. The restframe 
15 /xm luminosities of the BCGs in A1835 and Z3146 are an order of magnitude or more 
larger than those of the rest of the sample. A2390 is an intermediate case in that its 15 /xm 
luminosity is only slightly excessive (~ 5 o"), but the 70 /xm photometry clearly shows that 
it has a prominent far- infrared peak in its SED. 

These three infrared-brightest BCGs arc located in the three most X-ray luminous clus- 
ters in our sample, with soft X-ray (0.1-2.4 keV) luminosities of 16 — 21 x 10^^ erg s"^. These 
three clusters also have short gas cooling times of 0.6-1.9 Gyr based on the latest Chandra 
results (Bauer et al. 2005) as shown in Table 1. 

5.2. Radio Emission 

Ten of the eleven BCGs (i.e., except for AC114) have been observed in the radio (Ta- 
ble 2). Figure 5a plots the ratio of the flux densities at 24 fim and 20 cm parameterized as the 
^24 parameter (g24 = log(S'24^m/'S'20cm)) (Appleton ct al. 2004). The g24 values of star-forming 
galaxies fall in a well-defined range (0.84 ± 0.28) as a result of the radio-infrared luminosity 
correlation (Appleton et al. 2004), and therefore this parameter allows us to detect the exis- 
tence of excess radio emission likely produced by a radio AGN. Figure 5a shows that seven 
of the ten radio-observed BCGs are clearly out of the range for star-forming galaxies, having 
too much 20 cm radio flux density for their 24 //m flux density. The remaining three (those 
in A773, A2218, and A2219) are not detected in the radio with inconclusive lower hmits on 

?24- 

Based on the 20 cm luminosities, this excessive radio emission in the BCGs is hkely to 
originate from radio AGNs. Figure 5b shows that the 20 cm luminosities of these BCGs 
are mostly above 10^^'^ W Hz~^, which is usually the dividing line between AGNs and star- 
forming galaxies (e.g., Yun, Reddy, & Condon 2001). All seven BCGs outside the range of 
star-forming galaxies in Figure 5a fall in the radio AGN category in Figure 5b. The prevalence 
of radio AGNs in the BCGs means that there may be a significant AGN contribution to the 
observed infrared luminosities of these galaxies. 

The three BCGs with a low 20cm luminosity (< lO^''^-^ W Hz^^), A773, A2218, and 
A2219, arc all located in non-cooling fiow clusters with a radiative gas cooling time of longer 
than 10 Gyr as categorized by Allen & Fabian (1998) (Table 1). They have the longest 
cooling times (24.8, 30.9, and 29.1 Gyr) of the ten radio-observed clusters based on the 
latest estimates by Bauer et al. (2005) (Table 1). Such a trend is consistent with the earher 
finding that there exists a strong coupling between cluster cooling cores and radio emission 
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associated with the BCGs (e.g., Burns 1990), which may suggest that appreciable coohng 
of the intracluster gas in the cluster core tends to trigger a radio AGN in the BCG. The 
three clusters with no radio detection are also known to show bi/tri- modal dark matter 
distributions based on the Icnsing analysis by Smith et al. (2005). Such complex cluster 
mass distributions, together with the disturbed X-ray morphologies, may be related to the 
long gas cooling times. 

5.3. Properties of the Three Infrared-Brightest BCGs 

Here, we examine in more detail the properties of the three infrared-brightest BCGs in 
the sample, those in A1835, Z3146, and A2390. Various properties of these three BCGs are 
summarized in Table 3. 

Figiuc 6 shows the SEDs of the three BCGs. A number of photometric measurements in 
the literature are shown together with the Spitzer measurements. We modeled their behavior 
with a combination of the following three SED components: a giant elliptical (gE) SED, an 
infrared-luminous galaxy SED (Arp 220/M 82), and a power-law radio spectrum simulating 
the radio emission from a radio AGN. The gE SED was taken from the SED template library 
of the Hyper-z photometric redshift code (Bolzonella et al. 2000). For infrared-luminous 
galaxy SEDs, we used those of Arp 220 and M 82 taken from Silva et al. (1998). The SEDs 
of Silva et al. (1998) do not have the precision to reproduce individual spectral features (e.g., 
PAH features) , but provide accurate broad-band SEDs from ultraviolet to radio. 

The results are shown in Figure 6. As shown in panel (a), when the three components 
are combined, the gE SED dominates at Xrest ^ 5//m , the radio AGN SED dominates at 
Kest ^ 1mm, and the infrared-luminous galaxy SED dominates between 5 fj,m and 1 mm. 
The scalings of the gE SED and radio AGN SED were determined by the IRAC and radio 
measurements, respectively; the scaling of the infrared- luminous galaxy SED was set by the 
observed 24 fim flux density alone, and no effort was made to fit the far-infrared/submm 
measurements. As shown in the figure, the BCGs in A1835 and Z3146 have a SED steeply 
rising from 24 to 70/160 /xm, which requires an Arp 220-like SED to fit without violating the 
observed gE SED in the restframe visual/near- infrared. The BCG in A2390, on the other 
hand, has an SED significantly flatter between 24 and 70 //m, and in this case an M 82-like 
SED is more appropriate. 

Figure 6 shows that the three-component model works reasonably well to reproduce the 
observed SEDs of these three BCGs (the thin sohd lines). However, the flgure also shows 
that the use of the template SEDs from a small number of well-known galaxies (Arp 220, 
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M 82) is clearly too simplistic, and does not work in detail. For example, the far-infrared 
SED peak at ~ 100/xm seen with the BCGs in A1835 and Z3146 is clearly longward of that 
seen with Arp 220. 

To determine the infrared luminosities of these BCGs accurately, we also fitted the far- 
infrared/submillimeter SEDs with a two-component modified black-body model^ (the thick 
solid lines). These fits produced infrared luminosities of 7.3 x 10^^, 4.1 x 10^^, and 0.3 x 10^^ 
L0 for the BCGs in A1835, Z3146, and A2390, respectively. This also means that the BCGs 
in A1835 and Z3146 can be classified as LlRGs. Although Arp 220 is technically a ULIRG, its 
SED fits those of the BCGs in A1835 and Z3146 presumably because the infrared luminosity 
of Arp 220 (1.6 x IO^^Lq by Sanders et al. (2003)) is barely above the ULIRG threshold. The 
infrared luminosity of the BCG in A2390 is also comparable to that of M 82 (0.6 x IO^^Lq 
by Sanders et al. (2003)). 

The models in Figure 6 clearly show that the radio continuum flux is well above that 
expected from the Arp 220 SED. This is why these BCGs have small 524 values outside the 
range for star- forming galaxies (Figure 5a), indicative of a radio AGN. However, the radio 
SEDs of the BCGs in both A1835 and A2390 show that the slope of the radio spectrum is 
steep {a ~ 0.8, where Fi, oc i^~°), and as a result the infrared luminosity is dominated by 
thermal dust emission in both BCGs, as found in general for steep spectrum radio galaxies 
(Shi et al. 2005). In the case of the BCG in Z3146, we do not have enough radio measurements 
to constrain the radio spectrum, but it is already clear from the SED in Figure 6 that it also 
has a large far-infrared peak. 

Despite the presence of a radio AGN, all the BCGs show infrared SEDs typical of star- 
forming galaxies. In fact, none of the BCGs shows a power-law infrared SED typical of an 
AGN-dominated galaxy, and all of them show an excess of 8 fim emission, which is likely 
due to the 6.2 fim PAH feature in the passband at these redshifts. In the case of the BCG 
in Z3146, the existence of a strong 6.2 /xm feature has been conflrmed spcctroscopically with 
a Spitzer/IRS observation (E. Egami et al. 2006, in preparation), which indicates that star 
formation is a major contributor to the infrared luminosity. 

A number of arguments indicate that the infrared luminosities are produced by star 



■^Each modified black-body model was represented by the following analytical expression: F^, = B^{Td){l— 
^-{viva) -^vhere B^i^i) is the Planck function with a dust temperature T^, is the frequency below which 
the thermal dust emission becomes optically thin, and /3 is the dust emissivity power-law index. A good fit 
was obtained for the BCG in A1835 with = 40 & 110 K, ^/q = 3 THz (i.e., 100 /im), and /3 = 1.5. The same 
dust temperatures were adopted for the BCG in A2390, which does not have any far-infrared/submillimeter 
measurements to constrain the SED, while dust temperatures of 35 and 100 K were used for the BCG in 
Z3146. 
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formation. Prom the infrared luminosities themselves, we derive star formation rates (SFRs) 
of 125, 70, and 5 Mq yr"! for the BCGs in A1835, Z3146, and A2390, respectively, using the 
relation of Kennicutt (1998a). Similar star formation rates are derived from the extinction- 
corrected Ha luminosities measured by Crawford et al. (1999) with the resultant values of 40, 
47, and 5 M0 yr~^ (Table 3). The Ha-derived star formation rates are often lower than those 
derived from the infrared luminosities (e.g., a factor of three lower with the BCG in A1835) 
probably due to the finite slit width (l'/2-l'.'3) of the optical long-slit spectroscopy (Crawford 
et al. 1999) and/or internal dust extinction. The infrared luminosities from the restframe 
15 iim luminosities using the relation of Elbaz et al. (2002) for star-forming galaxies are 
9.5 X 10^^, 3.3 X 10^^, and 0.6 x 10^^ L©, similar to the values directly measured above. 

One major difference between the SEDs of these BCGs and those of infrared-luminous 
galaxies is that the former are much more luminous in the restframe visual and near-infrared. 
At Xrest < 5//m, the massive old stellar population in the BCGs completely overwhelms the 
light associated with the star-forming population. In fact, the total luminosity emitted in the 
visual/near- infrared from the massive old stellar population is comparable to that emitted 
by dust in the infrared. This explains why it has been difficult to detect the starbursting 
component in these BCGs. A similar object (i.e., a giant elliptical containing a starbursting 
population) has also been found in the field by Krausc et al. (2003). 

The two infrared-brightest BCGs, those in A1835 and Z3146, show another interesting 
property: exceptional strength of the CO (1-0) emission line for their infrared himinosities. 
Figure 7 shows a broad correlation that holds between the infrared luminosities and CO 
luminosities of LIRGs and ULlRGs. Plotted in the same figure, the BCGs in A 1835 and 
Z3146 hint at the possibility that they are overluminous in CO luminosity when compared 
with LIRGs and ULIRGs with similar infrared luminosities. This may mean that these two 
BCGs contain an unusually large amount of cold molecular gas. 

5.4. Implications for the Cluster Cooling Flows 

Figure 8 plots the 15 //.m monochromatic luminosities of the BCGs as a function of the 
radiative coohng time of the intracluster gas calculated by Bauer et al. (2005). Over the 
whole sample, there seems to be a trend that BCGs in clusters with shorter gas cooling times 
(and therefore larger X-ray luminosities in general) have larger infrared luminosities. Such a 
correlation would indicate that the infrared activity of a single galaxy at the cluster center 
is somehow connected to the property of a cluster core as a whole. 

Note, however, that a large X-ray luminosity does not always translate into a short 
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cooling time. For example, A2219 has a large X-ray luminosity but a long cooling time. In 
fact, A2219 is a peculiar cluster in that it highly deviates from the cluster X-ray luminosity- 
temperature {Lx vs. T) relation and mass-temperature (M vs. T) relation (Smith et al. 
2005). This cluster may have experienced a core-penetrating merger recently (Small et al. 
1995), and such a process might have boosted the X-ray temperature and luminosity (Smith 
et al. 2005). In any case, the fact that A2219 deviates from the overall trend in Figure 4 but 
not in Figure 8 suggests that the correlation between the 15 iim luminosity and gas cooling 
time is more fundamental than that between the 15 /xm luminosity and X-ray luminosity. 

In general, the increased infrared luminosities observed in local galaxies arc thought 
to be episodic. For example, transient galaxy-galaxy interactions are believed to boost the 
infrared luminosity by triggering vigorous star- formation and/or fueling an AGN. However, 
the process that triggered the increased infrared activity in the BCGs in A1835 , Z3146, 
and A2390 is difficult to explain in terms of similar processes since in this case the infrared 
luminosity of a BCG would likely increase independent of the properties of the cluster in 
which it resides. In other words, in this scenario there is nothing that would prevent a BCG in 
a low X-ray luminosity and long cooling time cluster from becoming infrared- luminous. This 
would introduce a significant scatter in the correlations of the BCG mid-infrared luminosity 
with the cluster X-ray luminosity (Figure 4) and with the cluster radiative gas cooling time 
(Figure 8). The fact that we do not see significantly deviant points from these trends 
suggests that the triggering mechanism is not a stochastic process but rather a process that 
is controlled by the properties of the cluster cores as a whole. 

An attractive possibility is to associate the observed infrared luminosity with the fate 
of the coofing cluster gas. In this case, only those BGGs in clusters with a short gas cooling 
time could become infrared-luminous because only these clusters could channel the gas into 
the BGGs and form stars. This hypothesis can be tested by comparing the infrared-derived 
SFRs with the X-ray-derived mass deposition rates. 

Table 4 lists various determinations of the mass deposition rates for A1835, Z3146, and 
A2390 based on the X-ray spectral analysis. Together with the previously pubfished values, 
we also list our own estimates based on the reanalysis of the Chandra data. We calculated the 
mass deposition rates at the cluster center within a radius of 13" (~ 50 kpc) from the BCGs 
(i.e., the size of the 24 /xm photometry beam) as well as for the entire clusters (r < 500 — 700 
kpc). The central mass deposition rates are more relevant for comparison with the star 
formation rates in the BCGs, and can be compared with other published C/iandra-derived 
mass deposition rates, which were also calculated within comparable radii (Table 4). We 
did not fist the mass deposition rates derived from the ROSAT and ASCA data because 
they are believed to be overestimated by almost an order of magnitude in many cases. The 
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large mass deposition rates resulting from the RO SAT and ASCA data, exceeding 1000 
yr~^ in extreme cases, are thought to originate from the existence of cool ambient gas in 
the cluster core, which, when spatially unresolved, tends to be included into the cooling- 
fiow calculation and inflates the mass deposition rates (Allen, Ettori, & Fabian 2001). The 
difference between the central (r < 50 kpc) and total (r < 500 — 700 kpc) mass deposition 
rates listed in Table 4 reflects the size of this effect. 

Table 4 shows that the various estimates of the mass depostion rates are not always 
consistent. In the case of A1835, for example, the value of 130-200 Mq jt^^ (r < 24 kpc) 
derived by Schmidt, Allen, & Fabian (2001) is signiflcantly larger than the other estimates. In 
fact, Voigt & Fabian (2004) point out that the existing Chandra data of A1835 at r < 100 kpc 
is statistically consistent with the absence of the cooling-flow component, and our reanalysis 
of the same data is in line with this result. In the case of A2390, the various estimates are 
broadly consistent with each other, especially if we account for the effect mentioned above 
that would inflate the mass deposition rate as we increase the region to extract the X-ray 
spectrum. In the case of Z3146, we do not have any other estimates to compare at this 
point, but this cluster provides the most convincing evidence for having a signiflcant mass 
deposition rate. The differences between the various estimates probably reflect the level of 
uncertainty inherent in such analyses, and some part of the differences may be attributed to 
the improvements in the Chandra/ ACIS calibration since the early years of the mission. 

The infrared-derived SFRs for the BCGs in A1835, Z3146, and A2390 arc 125, 70, and 5 
Mq yr^^, respectively, while the X-ray-dcrivcd mass deposition rates at the chistcr center are 
< 20, 290, and 110 Mq yr~^, respectively, based on our own estimates. This indicates that 
the correspondence between the infrared-derived SFRs and X-ray-derived mass deposition 
rates is currently suggestive at best. Z3146 is the only case where the two quantities are 
broadly consistent (within a factor of several). A2390 shows a large difference (a factor of 
20), which may be related to the strong radio AGN detected in the BCG of this cluster. In 
the case of A1835, the difference seems at lease a factor of six although we are in need of 
deeper Chandra data to better constrain the mass deposition rate and reconcile the conflicting 
estimates for the mass deposition rate. 

Given the complex nature of star formation processes and the limited accuracy with the 
determination of the physical quantities such as star formation rates and mass deposition 
rates, it would not be surprising that the infrared-derived SFRs and X-ray-derived mass 
deposition rates differ by a factor of several or more even if there is a causal connection 
between the two. For example, the well-known Schmidt law (i.e., the correlation between 
the surface densities of gas and star formation rate seen in external galaxies) has a 1 cr scatter 
of a factor of 2, and individual galaxies deviate by as much as a factor of 7 (Kennicutt 1998b). 
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Therefore, at least in Z3146 (and possibly in A1835), it is possible to interpret the observed 
SFRs as a result of the cooling cluster gas accreting onto the BCGs and triggering star 
formation (i.e., the cooling-fiow interpretation). Such a scenario is also consistent with the 
the exceptionally strong CO line emission observed with these two BCGs since cold molecular 
gas is seen to be accreting to the BCG in one of the low-redshift cooling-flow clusters (Salome 
& Combes 2004). 

6. Conclusions 

With the Spitzer/YRAC (3.6-8.0 /xm) and Spitzer/MIP^ (24-160 ^m) imaging data, we 
studied the infrared properties of the brightest cluster galaxies (BCGs) located in the cores 
of X-ray-luminous clusters. The main results are as follows: 

• The majority of BCGs are not infrared- luminous. Out of the 11 BCGs studied here, 
only three (those in A1835, Z3146, and A2390) have a strong far-infrared SED peak, 
and only two (those in A1835 and Z3146) can be classified as LIRGs (> 10^^ L0). 

• There seems to be a trend that clusters with larger X-ray luminosities and shorter 
gas cooling times harbor more infrared-luminous BCGs. The three infrared-brightest 
BCGs (those in A1835, Z3146, and A2390) were found in clusters with the largest 
X-ray luminosities and shortest radiative gas cooling times in the sample. 

• The infrared luminosities of the three infrared-brightest BCGs are likely to be powered 
by star formation judged from the infrared SED signatures. The infrared-derived SFRs 
are comparable to the extinction-corrected Ha-derived SFRs. 

• Radio AGNs are prevalent in the BCGs. Seven out of the ten radio-observed BCGs 
studied here contain a luminous radio source (> 10^^'^ W Hz^^ at 20 cm), and show a 
low 24 /xm/20 cm flux density ratio, both indicative of a radio AGN. 

• The infrared-derived SFR of the BCG and the X-ray-derived mass deposition rate in 
the parent cluster are broadly consistent in Z3146 (and possibly in A1835), suggesting 
the interpretation that the observed star formation is triggered by the cooling cluster 
gas accreting onto the BCGs except when there is a strong radio AGN as is the case 
in A2390. 

This work is based on observations made with the Spitzer Space Telescope, which is op- 
erated by the Jet Propulsion Laboratory, California Institute of Technology under a contract 
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with NASA (contract number #1407). Support for this work was provided by NASA through 
an award issued by JPL/Caltech (contract number #960785, #1255094). GPS acknowledges 
financial support from Caltech and a Royal Society University Research Fellowship. 

A. Reanalysis of the Chandra X-ray data 

The Chandra data were re-analyzed using CIAO 3.2 and the latest cahbration files 
available in CALDB 3.1. Each individual dataset was reprocessed and screened to remove 
periods of strong background flares in the standard manner. Matching background event 
files were created for each dataset from the standard ACIS blank-sky event files following 
the procedure described in Vikhlinin et al. (2005). 

For a given cluster, we have calculated the total gas coohng rate for the entire cluster 
as well as the cooling rate in the immediate surrounding of the BCG (within a radius of 13", 
which is the size of the 24 /xm photometry beam). The total cooling rates were calculated 
based on fits to the integrated spectrum inside a maximum outer radius. In practice, this 
radius was determined from the radial surface brightness and defined as the point at which 
the surface brightness profile drops below the background level. For A1835 and A2390, 
the individual datasets were combined to form a mosaic and the radial surface brightness 
determined for this composite. With this definition, the corresponding maximum radii are 
125" (~ 490 kpc), 125" (~ 460 kpc), and 160" (~ 700 kpc) for A1835, A2390, and ZW3146, 
respectively. 

To calculate the gas cooling rates, spectra in a given region were extracted in the 0.5-7.0 
keV energy range along with corresponding background spectra from the blank-sky event 
files mentioned previously. Counts weighted detector response (RMFs) and effective area 
(ARFs) files were created for the extraction regions using the CIAO tools mkacisrmf and 
mkwarf . The resulting spectra were then compared to various standard spectral models using 
the XSPEC model library available in the ISIS spectral fitting package (Houck & Denicola 
2000). The fitted spectral model was defined to consist of foreground Galactic absorption, a 
single temperature thermal plasma (MEKAL) due to emission from the hot, outer cluster, 
and an isobaric coohng flow component (MKCFLOW). Fits were performed with the Galactic 
column flxed to the standard value as well as allowing it to varying. For the two clusters 
with multiple observations, individual spectra were extracted for each dataset and the model 
was flt to the component spectra simultaneously. For A1835 and ZW3146, the flts with and 
without varying Galactic columns gave comparable values and the cooling rates quoted 
in Table 4 are for the flxed column flts. Fits for A2390 showed a marked improvement in 
if the Galactic column was increased by ~30% and the result from this flt is quoted in 
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Table 4. 
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Table 1. Target Clusters and Their X-ray Properties 



Name'' RA'' DEC'' X-ray Luminosity Gas Cooling Time Typei 











ROSAT (NORAS)<i 
(lO^^'erg s-i) 


ROSAT (BCS)" 
(10*4 erg s-i) 


ASCA* 
(lO^^erg s-i) 


Einstein (EMSS)^ 
(10^4 erg s-i) 


ROSAT^ 
(Gyr) 


Chandra' 

(Gyr) 




A773 


09 17 53.4 


+51 43 39 


0.217 


7.66 


8.27 


9.36 




9.1 


24.8 


NCF 


Z3146 


10 23 39.6 


+04 11 12 


0.291 


20.18 


17.86 


24.89 




1.6 


0.6 


CF 


MS1358.1+6245 


13 59 50.6 


+62 31 05 


0.327 


6.14 




7.51 


7.38 


2.7 




CF 


A1835 


14 01 02.0 


+02 52 45 


0.252 


21.24 


25.14 


29.68 




1.4 


0.6 


CF 


MS1455.0+2232 


14 57 15.1 


+22 20 35 


0.258 


9.00 


8.65 


10.89 


10.51 


1.0 




CF 


A2218 


16 35 49.3 


+66 12 45 


0.175 


4.96 


5.66 


6.57 




9.5 


30.9 


NCF 


A2219 


16 40 19.6 


+46 42 43 


0.228 


15.54 


13.03 


21.84 




8.4 


29.1 


NCF 


A2261 


17 22 27.2 


+32 07 58 


0.224 


13.11 


11.57 


15.21 




2.8 


3.0 


CF 


MS2137.3-2353 


21 40 15.1 


-23 39 39 


0.313 






11.68 


10.73 


1.1 




CF 


A2390 


21 53 36.7 


+17 41 45 


0.233 


16.13 




26.49 




3.9 


1.9 


CF 


AC 114 


22 58 48.4 


-34 48 09 


0.312 






12.01 




16.2 




NCF 

1 



''The prefix "A" denotes the Abell clusters (Abell, Corwin, & Olowin 1989) while the prefix "MS" denotes the Einstein Extended Medium Sensitivity Survey clusters (Gioix? 
& Luppino 1994). The Zwicky cluster Z3146 is often designated as ZwC11021. 0+0426 (NED convention) while the cluster AC114 is also known as the Abell supplementarjy 
southern cluster AS 1077 (Abell, Corwin, & Olowin 1989). 

''These coordinates indicate the positions of the BCGs. The BCG position in MS2137 was measured in the Palomar Digital Sky Survey image while other positions were 
taken from Small et al. (2001) (A2218), Stanford et al. (2002) (MS1358 and AC114), and Crawford et al. (1999) (for the rest), respectively 

''AH the redshifts were taken from Allen (2000). 

^0. 1-2.4 keV luminosity from NORAS survey (Bohringer et al. 2000) 
«0. 1-2.4 keV luminosity from BCS survey (EbeUng et al. 1998) 
*2-10 keV luminosity from ASCA (Allen 2000) 

^0.3-3.5 keV luminosity from the Einstein Medium Sensitivity Survey (Gioia & Luppino 1994) 
''Radiative gas cooling times at the cluster center from Allen (2000) based on the ROSAT observations. 
'Radiative gas cooling times at the cluster center from Bauer et al. (2005) based on the Chandra observations. 



^Cooling flow (CF)/Non-cooling flow (NCF) classifications by Allen & Fabian (1998) (CF: 90%-conficence upper limit on tcool < 10 Gyr). The Chandra cooling times by 



Bauer et al. (2005) do not affect the classification at least for these clusters. 



Note. — All the listed values correspond to the cosmological parameters {CIm, f^Ai ^^0 (km s ^ Mpc ^)) of (0.3, 0.7, 70). When converted from a different choice of the 
cosmological parameters, the X-ray luiiiiriositics were scaled with the inverse square of the luiiiinosity distance. Since the gas temijcraturc is independent of the choice of the 
cosmological parameters, the gas cooling time through bremsstrahlung emission scales with n^^ oc ^ w^^""® "ei ^Pt are the electron density, proper volume, and 
luminosity of the X-ray emitting gas. 



Table 2. Spitzer/VLA Flux Density Measurements of the BCGs 



IN dlilt; 


O.U fjiLlL 


4.0 

Cm Iv^ 


O.O fJjUl 
Cm Iv^ 


Q n /;m 

Cm Tir"! 


OA » im 
^111 

Cl1 l-iA 


7n /im 

Cm Iv^ 


Ififl iim 
XUU ^111 

Cm l-iA 


1 A 

ill J y lei 


A773 


1.42 


1.10 


0.62 


0.42 


98 ±14 






< 0.14*^ 1 


Z3146 


0.89 


0.80 


0.54 


1.33 


4099 ± 410 


68 ±14 


157 ± 35 


2.04 2 


MS1358 


0.69 


0.56 


0.32 


0.24 


90± 13 






2.61 2 


A1835 


2.17 


2.03 


1.42 


5.14 


17236 ± 1724 


185 ± 37 


305 ± 63 


31.25 2 


MS1455 


1.27 


1.00 


0.59 


0.50 


345 ± 37 






5.62 2 


A2218 


1.18 


0.86 


0.52 


0.31 


84± 13 






<0.14'* 1 


A2219 


1.02 


0.79 


0.42 


0.28 


125 ± 15 






< 0.14^ 3 


A2261 


2.15 


1.66 


0.90 


0.65 


347 ± 36 






3.40 2 


MS2137 


0.88 


0.68 


0.42 


0.28 


294 ±31 






3.8 4 


A2390 


1.08 


0.87 


0.53 


0.67 


1075 ± 108 


8±3 




200 3 


AC 114 


0.95 


0.79 


0.48 


0.30 


160 ± 18 









^^3 a upper limit {a 45/iJy) 



References. — (1) Morrison et al. (2003); (2) FIRST survey (Becker, White, & Helfand 1995; White et 
al. 1997). The latest catalog (2003 Apr 11 version) is available from http://sundog.stsci.edu; (3) Rizza et al. 
(2003); (4) NVSS survey (Condon et al. 1998) 

Note. — All the BCGs are strongly detected in the IRAC bands, and therefore the IRAC photometric 
uncertainties arc dominated by those of the absolute calibration. For the MIPS bands, the 1 a error bars 
including both the measurement and calibration uncertainties are explicitly shown because the contribution 
from the former is not negligible for weak detections. 
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Table 3. Properties of the Three Infrared- Brightest BCGs 





A1835 


Z3146 


A2390 


LiR ^ (IO^Lq) 


7.3 ± 1.5 


4.1 ±0.8 


0.3 ±0.1 


SFR (IR)*^ (Mq yr-i) 


125 ±25 


70± 14 


5± 1 


L(Ha)'= (lO'^o erg s^^) 


511 ±105 


594 ± 59 


67± 18 


SFR (Ha)d (M© yr"!) 


40 ±8 


47 ±5 


5±1 


ScoAl/" (J km s-i) 


7.7± 1.3 


5.2 ± 1.2 




L'co^ (101° K km s-i pc2) 


2.4 ±0.4 


2.2 ±0.5 





"^Infrared luminosity calculated from the two-component modified 
black-body SED models shown in Figure 6. 

''Star formation rate derived as 4.5 x 10~^^i/i{ (erg (Kennicutt 
1998a). 

'^Extinction-corrected Ha luminosity from Crawford et al. (1999) 
rescaled with = 0.3, = 0.7, and Hq = 70 km Mpc"'. 

'I Star formation rate derived as 7.9 x 10~^^ Z/(iIc«) (erg s~^) (Kennicutt 
1998a). 

"Velocity-integrated line flux measured by Edge & Frayer (2003) 

^CO(l-O) line luminosity calculated as i'^o = 3.25 x 
lQ'^Sco^Vv^j^^D'l{l + z)-^, where Sco^V is the velocity-integrated 
hne flux in Jy km s~^, Uobs is observed frequency in GHz, and Dl is 
the luminosity distance in Mpc (Solomon et al. 1997). 

Note. — For Lj/j and SFR(IR) (i.e., the first two rows), we adopted 
a conservative uncertainty estimate of 20%, which corresponds to the 
absolute calibration uncertainties at 70 and 160 /xm. 
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Table 4. Mass Deposition Rates for the Three X-ray-Brightest Clusters 



Cluster 






Chandra 






XMM 




SAFOl^ 
(<24 kpc) 


AEFOl^ 
(<40 kpc) 


This work*' 
(<50 kpc) 


VF04= 
(<100 kpc) 


This work<i 
(<500-700 kpc) 


POl 

(< 120 kpc) 


A1835 
Z3146 
A2390 


130-200 


130-190 


< 20 


339tt?8 


100+r 

670tl?° 
42011° 


< 130 



'^We adopted the mass deposition rate of 200-300 Mq yr ^ at r < 30 kpc for A1835 and at r < 50 kpc 
for A2390 {Qm = 1, Oa = 0, and Hq = 50 km s'^ Mpc-^). 

''The centeral mass deposition rates derived for the r < 13" region, which is the size of the 24 fim 
photometry beam. This radius corresponds to 51, 57, and 48 kpc at the distances of A1835, Z3146, and 
A2390, respectively. 

■^The deposition rates were calculated inside the region where the gas cooling time is less then 5 Gyr. 
This roughly corresponds to 100 kpc for A1835 and A2390. 

■^The total mass deposition rates calculated within the maximum radii of 490, 460, and 700 kpc at the 
distances of A1835, A2390, and ZW3146, respectrively. 

Note. — All the numbers are in the unit of Mq yr~^, and the error bars indicate the 90% confidence 
limits. The numbers in the parentheses indicate the radius of the region inside which the mass deposition 
rate was calculated. The mass deposition rate scales with the inverse square of the luminosity distance, 
and all the published values have been rescaled with Qm = 0.3, I^a = 0.7, and Hq = 70 km s~^ Mpc~^. 
The XMM-derived mass deposition rate for A1835 by Peterson et al. (2001) was rescaled with a luminosity 
distance of 1272 Mpc. 

References. — SAFOl: Schmidt, Allen, & Fabian (2001); AEFOl: Allen, Ettori, & Fabian (2001); VF04: 
Voigt & Fabian (2004); POl: Peterson et al. (2001) 



-27- 



A1835 




Fig. 1. — Spitzer images of the cluster cores harboring the three infrared-brightest BCGs. 
From left to right are the 3.6, 8.0, 24, 70, and 160 fim images, and the clusters are ordered 
vertically in decreasing 24 fim flux densities. North is up and east is left. Each image is 1' 
on a side, and the pixel scales are 0'.'6 pixel^^ at 3.6 and 8.0 /xm, 1'.'25 pixel"^ at 24 /xm, 4'.'92 
pixel^^ at 70 /xm, and 8" pixel^^ at 160 /xm. 
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MS1455 




Fig. 2. — Spitzer images of the cluster cores for the rest of the sample. From left to right are 
the 3.6, 8, and 24 /im images, and the clusters are ordered vertically in decreasing 24 ^m 
flux densities. North is up and east is left. Each image is 1' on a side, and the pixel scales 
are same as those in Figure 1. 
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Fig. 2. — continued. 
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Fig. 3. — Spitzer SEDs of the BCGs. All the SEDs were scaled such that the flux densities 
coincide around 4.5 iim. The SED of the BCG in A1835 is denoted as a dashed line to make 
it distinguishable from that of Z3146. 
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Fig. 4. — The 15 /xm monochromatic luminosities (z/Lj^) of the BCGs are plotted against the 
cluster soft X-ray (0.1-2.4 keV) luminosities from the NORAS survey (Bohringer et al. 2000). 
For the two southern clusters MS2137 and AC114 (unfilled circles), we used the Einstein 
0.3-3.5 keV (Gioia & Luppino 1994) and ASCA 2-10 keV (Allen 2000) X-ray luminosities, 
respectively. The latter was divided by a factor of 1.25 to account for the systematic offset 
between the ROSAT soft X-ray luminosities and ASCA hard X-ray lumninosities. The 
horizontal dotted lines show the 1 a range of the 15 //m monochromatic luminosities observed 
for local eUiptical galaxies (Ferrari et al. 2002), which is i/Lj,[15//m] = (1.7±1) x lO^L© (Note: 
these 15 //m luminosities are not those tabulated by Ferrari et al. (2002), which do not contain 
the stellar contribution, but those directly calculated from their flux density measurements.) 
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Fig. 5. — (a) The ^24 parameter (g24 = log(5'24^m/'S'20cm)) of the BCGs plotted against their 
15 fim monochromatic luminosities. The shaded horizontal band shows the 1 a range of 
the ^24 parameter observed for the 24 /xm-detected star-forming galaxies (Appleton et al. 
2004); (b) The 1.4 GHz radio luminosities {L,y) of the BCGs plotted against their 15 iim 
monochromatic luminosities. The radio luminosities were calculated as L = ^t^D\Siaghz 
(i.e., applying a K correction assuming a radio continuum shape of oc 
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Fig. 6. — SEDs of the three infrared-brightest BCGs. The sohd circles indicate the Spitzer 
measurements while the open circles indicate the measurements in the literature (the error 
bars are comparable or smaller than the symbols). The dotted line shows the SEDs of 
Arp 220 (a and b) and M 82 (c) (Silva et al. 1998). The thin solid line shows the three- 
component SED model while the thick sohd line shows the two-component modified black- 
body SED model to fit the observed far-infrared/submm SED (see the text for details). 
The flux density measurements of the BCGs in A1835 and A2390 were taken from Edge et 
al. (1999). The visual/near-infrared measurements were omitted because of the unknown 
aperture corrections, but they show a sharp drop toward shorter wavelengths typical of 
the SED of a giant elliptical. For A2390, we also excluded the low signal-to-noisc IRAS 
measurements. The measurements for Z3146 were taken from Chapman et al. (2002) (850 
lim) and the FIRST survey (20 cm). 
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Fig. 7. — The infrared luminosity and CO(l-O) line luminosity of the BCGs in A1835 and 
Z3146 (sohd circles) are compared with those of LIRGs and ULIRGs. See the notes of Table 3 
for the definition of the CO line luminosity, Lqq. The LIRG/ULIRG points shown are from 
Sanders, Scoville, & Soifer (1991) (triangles), Solomon et al. (1997) (crosses), and Gao & 
Solomon (2004) (open circles). The infrared luminosity and CO line luminosity of the two 
BCGs listed in Table 3 were rescaled to the values with Hq = 75 km s~^ Mpc~^ and go — 0-5 
so that the direct comparison with the LIRG/ULIRG sample is possible. 
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Fig. 8. — The 15 /xm monochromatic luminosities of the BCGs are plotted against the 
radiative coohng times of cluster gas derived from Chandra observations by Bauer et al. 
(2005) (solid circles) and from 05^ T observations by Allen (2000) (open circles). The 
horizontal dotted lines show the 1 a range of the 15 //m monochromatic luminosities observed 
for local ellipticals shown in Figure 4. 



